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ABSTRACT 



A two-diraonsional omega equation ia derived by 
combination of the vorticity and thermodynamic equationa. 

The deaired omega ia then taken to be the logarithmic average 
in the 1000-700 mb layer. A diabatic term, after Laevaatu, 
for oceanic areaa only la included to deduce the empirical 
temperature and vapor-preaaure changea aaaociated with 
aenalble and latent heating in the maritime layera. Over 
both continental and oceanic areaa a frictional vorticity 
aink ia included in order that exceaaivo energy cannot be 
generated over the ocean. Among other novel featurea ia the 
uae of the Holl atatic-atability parameter which afforda 
vertical conaiatonoy wit!i analysea prepared by Fleet 
Numerical Weather Facility. 
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1. Introduction 

Little practical us© has been made of vertical motions, 
in short-range forecasting due to the complicated 
and lengthy computations required. Vertical motion computa- 
tions are usually the by-product of multilevel baroclinic 
models and single values must be extracted only after going 
through the entire three-dimensional procedure. In this 
paper a two-dimensional vortical motion equation is developed 
using vertically integrated parameters with an assumed vertical 
motion profile. The layer-mean omega is ©specially useful 
for short-range thickness forecasts and the 1000-700 mb layer 
has boon selected for study. Derivation of the CJ -equation 
is similar to that by Thompson [X8] except a diabatic heating 
term providing a mechanism for development over oceanic 
areas has been included. Duo to the empiricisms employed 
the non-olliptic conditions mentioned by Pedersen [X6] do 
not arise, 

2, The Basic Omega Equation 

The diagnostic omega equation is derived by standard 
means from a vorticlty equation and a thermodynamic equation 
which retains the diabatic term. Equation (1) is the time 
derivative of the first law of thermodynamics in (x, y, p, t) 
coordinates. 



iL4. a 

it 0 Cp 



( 1 ) 



Substituting, T • which comss from the hydrostatic 

assumption and the ideal gas lawj then dividing by, -g/R, 
yields equation (2), 






( 2 ) 



Operating on equation (2) with the Laplacian gives the final 
form of the thermodynamic equation as shown by equation (3). 



I %Cf 

Here, is the Roll stability parameter to be 



i 

diiuf 



l0dp 

further discussed in section (I[), 



The vorticity equation in pressure coordinates is shown 
by equation (4). 






(4) 



As presented by Thompson and numerous other writers 

the last term of the left side of (4) is approximately equal 

to that on the right side, and the two terms are henceforth 

deleted. Then, a useful form of the vorticity equation is 

arrived at (equation 5) by making the geostrophic assumption 

for vorticity and for velocity, - _2_Va , and by 

° ' f ’ f 

taking the logarithmic pressure derivative. 
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(5) from (3), 



is the absolute vorticity. Next subtract 
and the result is the omega equation: 




3. Vertical Distribution of Pressure and Omega 

Instead of the distribution of pressure indexes normally 
considered, that employed here is given by, 

where y) is an arbitrary pressure level. The 1000-500 mb 
layer is divided into sub-layers as shown by Figure (1), 

The resultant pressures, while essentially logarithmic, bear 
values nearly identical to those of the mandatory levels. 



4 — = P^2 ' = 500.0 mb 

3 p^ = Po2 ^ = 594.6 mb 

2 — - P^2~^ - 707.1 mb 

1 p, ^ p ^2 H = 840.9 mb 

0 P^ = Po2‘^ = 1000.0 mb 



Figure 1. Pressure Distribution in the Vertical 
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The values of 60 = dp/cti:, the "vertical velocity in pres- 
sure coordinates , ,ls assumed to vary parabolically in the 
vertical with a profile defined by equation (7)^ which 
extends nearly to the ^00 mb level. 






(7) 



For boundary conditions it is assumed that^=Oat p=p^ 

(the kinematic boundary condition at the assumed surface of 
the earth, p ). Also for convenience it is assumed that at p 

T 

i.e, that the 500-mb divergence is zero. Then, 
upon differentiation of equation (7) with the ^OO-mb diver- 
gence taken to be zero, equation (8) results. It should be 
noted that if the value is assumed zero at p=pj 

the value of B is altered by only 5^, 

Therefore, B=1,3863A, and the omega profile in terms of A is 
given by equation (9). 






or 



CO = A 



/AL\ + 

pj ^ V /?. 




(9) 



Since, =o, the ratio between 60,, and 60^, defines the 

omega profile (equation 10) for the 0-2 layer (see Figure 2). 
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cov » „ I, ( 10) 

<^l -^A^2('-;LXt2.4y.3863) 




Figure 2. The Omega Profile 



For the pressure-range considerod (see .Pigu.re 2), the layer- 
logarithmic -mean, C<J =0.950dOpSO that , will be used for<U . 
Note finally, the profile assumed here applies only to the 
large-scale adiabatic, frictionloss component of the vertical 
velocity (that is terrain irregularities are considered 
absent ) . 

It. Vertical Integration of the Omega Equation 

The Holl stability parameter [ii] is used since it is one 
of the vertical— consistency requirements used by FNIIP (Fleet 
Numerical Weather Facility) and this study uses FNWF processed 
data. The Holl parameter is a modification of the standard 
stability parameter used here, as shown by equations (11), (12), 
and (13). 
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(j^ _ Rpcr - R p T _ )^p /)RL _ 

e^>p ' \^fp 



( 11 ) 



Letting, RT = - ^ 



gives, 



C-H 4-R^ 



( 12 ) 



Equation (12) may now be finite differenced over the 0-2 layer. 



Oh > 

A.h 

Equation (13) represents the Holl stability parameter in 
finite-difference form. In effect, O h , by (13) gives the 
'logarithmic- pressure ■'aVeragb over the 0-2 layer. This 
use of C7^ is compatible with the assumed logarithmic pressure 
distribution and will be used throughout the analysis which 
follows. 






(13) 



The omega equation (6) will now be integrated term by 
terra over the isobaric layer (Po, P 2 ). Prom (6), one may 
rewrite the QJ -equation as: 



(A) (3) 

W (<Th oj) . f ))p"^ 





whose parts (A),,,.,, (S) will now be discussed individually. 
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(a) The first terra is approximated hy substituting 
layer raean values; , CQ, ^ CO 

(B) With the raatheraatical relation, ^ ^ ^ ^CO . 

b(ApT 

the second terra becomes, . \ ' 

whose center-f inite-dif f erenced form over the 0-2 layer 
becomes : 



^ f Ck^ 2 .~ 




(lit) 



Then, using the specified vertical motion profile the rela- 
tion^ 2 .= 1.71l43^i> is substituted. However, it is desired 
to introduce lower boundary effects due to friction and 
terrain by letting 0^d9.CO\.o , the ^lo terra then being pre- 
determined and placed on the forcing function side of the CO- 
equation. Expression (l5) below is the final form of (ll+) 
including effects of the lower boundary. 



Term (B) i | 

A Vpi / 

The lower boundary, CJto > will be discussed in section (5)» 



(C) Term (C) in the geostrophic-diagnostic model repre- 
sents geostrophic advection of absolute vorticity. Averaging 
in the vertical with respect to the logarithm of pressure leads 
to equation (16); 
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( 16 ) 




Th® derivation of equation (16) makes use of the mean value 
theorem, so that p = is represented by p^ and may be re- 
moved from within the integral. Equation (16) then becomes. 




(D) Term (D) involves the geostrophic advection of 
thickness within the logarthmic pressure integral. Its 
value is well approximated by equation (18): 




tbe right side of which gives equation (19). 

h£v''J(zi ,Aa) * ei-2o) (19) 

Equation (19) shows that the geostrophic advecting wind mnj’’ bo 
taken arbitrarily as the level (1) wind. 

(E) The diabatic heating term will be vertically 
integrated in section (6), 

The integrated form of the resulting omega equation, subject 
to the lower and upper boundary conditions (atp^andp<^) 
becomes. 
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l) 

-ft), 



H = 



60 



_L -^y)i &2g57<^, 

d«PV(^ 




4. £i2!v<ir(£,, - p 

T/LPVfr 



( 20 ) 



Equation (20) is a HeltnholtZ“typ© equation in the variable, 
^40^ , This grouping, , precludes making the usual 

simplifying approximation of most three-dimensional models, 
(see for example, Haltiner et al j^9 

V*(v w) = <rv*o) + cy\v V + z. W(r- w<;o i vntIo 



The final solution of equation (20) is to be divided by 
leaving ^ 1 = ^ , where ^ is the mean vertical motion which 
approximates the 8^0 mb vertical motion. 



5 . The Lower Boundary Condition 

Vertical motion at the lower boundary is the sum of 

a terrain and a friction term, 60 )/.© a • 

Here, COy , the terrain-effected vertical motion may be de- 
scribed by equation (21) from Berkofski and Bertoni j 2 j , 



U>r 






\72x 



( 21 ) 



Terrain height, is a smoothed field used operationally 
by PNWP, The wind velocity at terrain level, IVt- , will be 
arrived at by an objective procedure to be described below. 
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Th© friction term^p, as shown by equation (22) is a 
simplified version of Cressman's formula used by Haltiner 
et al 9 J . 

CO^ s C2 x>\/t^T (22) 

^I> is the geostrophic drag coefficient derived by Cressman D]« 
the field of which is in regular use at PNV/P, The field of 
is a set of constants, one for each (■•j) grid -point in th© 
Northern Hemisphere. 



Since the 7“ subscripts infer application at terrain height, 
values are used that approximate the particular terrain heights 
involved. This is done by dividing th© 0-2 layer into three 
terrain-height dependent cases. 



1 _7P0_ m 



I _J00_m 



-2200 m 

P. 

-^00 m 



Case 3 s 2200 m ^ Z-r 
use Pi. values 
Case 2: 900 m^Z'7^2200 m 

use P| values 



" Pa Case 1; sf c36’2t^ 900 m 

us© P© values 



Figure 3. Three Cases of the Lower Boundary 



Division of the 0-2 layer as shown by figure (3) is 
dependent upon the gradient level above terrain, since is 
approximated by some geostrophic wind throughout. As a basis 
for th© selection of the particular geostrophic wind, , note 
that the gradient level in a neutral atmosphere occurs at 



approximately TOO m above t er r ain level for a wide range of 
surface roughness^ as shown by Blackadar W- Therefore, 
vertical velocity at the lower boundary is case (1^ if the 
terrain height is 900 m or lower. This values corresponds to 
a "standard atmosphere" height which is TOOta or more below 
the (1) level, and (0) level parameters are then used to 
compute the lower boundary. Case (2) Is to bo used if the 
terrain height lies within a gradient-level range of pj 
but not . Case ( 3 ) occurs for all terrain heights within 
TOO m of , 

Equation (23) is the terrain-induced vertical velocity at 
the lower boundary after application of the geostrophic 
assumption. 




Here, the K1 subscript refers to the levels 0, 1, and 2, 
Densities are standard atmosphere values for the respective 
levels , 

6, The Diabatic Heating Term 

Diabatlc heating in this model results from an exchange 
of heat between the sea surface and the atmosphere. The ex- 
change is divided into three cases; neutral, convective, and 
inversion, the division being dependent upon the sea-air 

temperature difference, T - T , 

w a 

The neutral case (0* C ^ T^^-T^< 3°C ) results in no net 
transport of heat across the air- ocean interface and the 
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diabatic heating term is arbitrarily taken as zero when this 
condition prevails. 



The convective case ( results in a net heating 

of the atmosphere. The convective atmosphere is assumed to 
consist of a dry adiabatic lapse rate to the condensation 
level and a moist adiabatic lapse rate above. This particular 
model of the convective case is discussed by Burke [s] in his 
paper on the transformation of cP to mP air. Heat gain by the 
atmosphere is manifested by a layer-mean temperature gain 
induced by a one-hour trajectory of surface air over a warmer 
sea-surface, and the layer temperature change is shown by 
figure (1|). 




p. 



P, 



P, 



Figure 1|, Model of the Convective-Atmosphere Modification 



Heating is composed of: (1) sensible heat from the sur- 

face to the condensation level and (2) latent heating above 
this level. This result follows by an analysis of the Burke 
model and its implications reguardlng heat transport, Martin 
M has shown that the large-scale transports of heat have 
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already been included in those terras of equation (6) in 
which Q does not specifically appear. Thus h/Cp , for use 
in equation (6) is given by the local temperature change, 

dia ^ j (2ii) 

\ /kt /sCAJSiacS V /cATeWT 

Equation (2^) is the diabatic heating terra (E.) r o t= x> o " 



to Ira o'uo.tto"' (6). 



(2^) 



Integrating logarithmically in the vertical yields the 

final form of the diabatic term as shown by equation (26), 
/^C. /'H 



AvFVa y?2_ ^?c. 



^ ^ 4 . (26) 

' A-t A-i 

In each of equations 2l|, 25# and 26, the operator-^i- represents 



At 



the finite-difference version of the local time derivltive, 
while, , represents a convective-condensation level 
pressure. 



According to the model depicted in figure (i|), the modi- 
fication is largely determined by the static stability which 
in turn largely depends on<fff , As an assumed infinite source 
of heat, the ocean will first modify a thin surface layer of 
air and then, due to convective activity, will distribute the 
modification (temperature change) throughout the 0-2 layer. 

An overestimation of heat exchange should be expected with 
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this particular process since convection is not instantaneous 
This overestiraation may be partially balanced by heat lost 
from the top of the 0-2 layer. However, the non-inclusion of 
a moisture-continuity equation may be more serious in regions 
of upper ridges. 

The inversion ca se (T^^ - 0) or cooling cas^„ whore 

heat is lost from the atmosphere to the ocean is defined 
by an empiricism from Laevastu 

” 3.0 V/0 {TIo— ~Ta} gm-cal/cm^-24hrs (27) 



Equation (2?) represents the heat flux across the 10-meter 
level. This is an empiricism requiring the 10-meter wind, 
V)o, which will be approximated by making a frictional 
correction to the 1000 mb geostrophic wind as shown by 
equation (28): 



= o.iA. 
■f 






(28) 



The factor 0.5 is due to surface frictional effects and will 
be discussed in section (7), 



Before integrating equation (27) over the 0-2 layer 
the initial distribution within the layer must be determined. 
Normally cooling effects due to inversion conditions are con- 
fined to a surface layer only several hundred meters thick, 
however, since the layer mean cooling effect is desired over 
the entire 0-2 layer it will be assumed that cooling is dis- 
tributed throughout the 0-2 layer. Thus, dividing equation 

24 



(27) by the mass of a column of air extending from the surface 
pressure tc with total mass,>v» = (p^ - Pi)/g» gives the 
mean temperature change in a vertical column extending from 

P^ to Pj_ as. 



Q - ^ 



(29) 



Finally wo have, 

cp GpA.%4 CBrfi)P 



(30) 



which is the final result for the inversion case assuming 
negliblo heat release from fog-droplet condensation. 



7. State-Change Parameters for the Convective Case 

Air-temperature change and wet-bulb temperature change 
equations must be derived for use in equation (26). The 
state-change equations (31) and (32) were Initially developed 
by Mosby Amotjjl'j, and more recently by Boyumj^l. 

However, in a PF/JF study by Carstensen and Laevastujjs'j the 
following best-fit equations were found to give hourly changes 
with a high degree of skill: 



=. a/'Z.(^'77u-'7A)-OJO -Wio^VTa 

At 


°CAr 


(31) 


- oj5^ew-eA)- o.ia - 


Wb/hr 


(32) 



Equation (31) is substituted directly into the sensible- 
heat term of equation (26) while the wet-bulb temperature 
change contributes directly to the latent-heat term. In 
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computing ^2^8 it is necessary to include the small increment 

At 

of temperature found by proceeding along a mixing ratio isopleth 
from to then down a moist adiabat to Pq, (see figure 

5 ). Equation ( 33 ) shows this relationship. 






Here A?pc is the one-hour height change of the condensation 
heights^ determined by the height of intersection of the dry 
adiabat and the mixing ratio isopleth based upon the tempera- ' 
ture - dew-point spread at the surface, 

-2l*ci si(~fAT.-Th2.\ (3h) 




Substituting equation ( 34 ) into (33) gives the local rate of 
change of the wet-bulb temperature as a function of dew-noint 
and ai r-tornnei^n-hur o chonfro'i j , 

^1k>b = X _ Ah>\ (35) 

At At \h-^4 J\A^t At ) 
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The dew-point temperature change is found by taking 
the time differential of equation (36) below [xo], 

(36, 

Here, e^ is the saturated vapor pressure and e^ the observed 
value, both at 10-meters, The finite-difference form of the 
time differential of dew-point resulting from equation (36) 
then follows; 



At At 



The vapor pressure change, , is given by equation (32), 



The time rate of change of saturated vapor pressure. 



AS 






is calculated using the Cl.ausius-' Olape^^ron equation: 



- o,^Z2^L<e^ ATa 
A-t 



(38) 



= 4Jo9<fxj> 






“7a . 



(39) 



Equation (39) is an integrated form of the Clausius - Clapjrron 



equation. 



There are still some unsolved parameters in the system. 
The moist adiabatic lapse rate, , cannot be considered a 
constant with respect to pressure even in the limited opera- 
tional range of this study. At the condensation level, , 
the moist adiabat is calculated from the following rio] : 
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= UAi — ^iL — 






iho) 



All subscripts G refer to the condensation level. 

Condensation parameters of temperature f and pressure, 

; 'ir.; >‘ ‘'1 

Pj, , are readily calculable from equation (4D and (i).2) after 
Ed son L8]: 



7i=i; _ 



?c = R> 




("Ta - Ttj 






( 1 * 1 ) 



( 1 * 2 ) 



For equation (1^0), the vapor pressure at the condensa- 
tion level, e^, must be calculated. Equation (1|3)» the inte- 
grated form of the Clausius -Glapeyron equation may be used 

'fci' 

since air at the condensation level is saturated, and and 
are known from (1|1) and ( 42 ), 






273. /(^ 



X. 



(43) 



In like manner the vapor pressure of the sea surface, e^, 
(from equation 32) is calculated since the sea surface is 
saturated (equation 44). Also, a correction factor of 
0,98 is required for salinity effects. 



euo = (o.^i)(^io?)sKf> 



6^1 f-o 5V/P^O 

Z93./4 ~Tu) 



(44) 



.JO 



Finally, the advecting wind at 10-meters, ^le> , must be 
derived for use in equations (31) and (32)* The geostrophic 
wind at the surface is taken as equal to the 1000 mb geo- 
strophic wind. However, due to frictional interaction with 
the ocean siirface, the geostrophic wind must be reduced in 
magnitude and rotated to the left, to yield the advecting 
frictional wind,lV|^, For the convective case the observed 
cross-isobar angle at mid-latitudes is approximately 15^ 

Reduction of the 10-raeter wind speed as expressed by 
the ratio, V/Vgo , where is the 10-meter geostrophic wind 
speed, may be obtained from a study by Lettau L«l- By 
assuming a roughness length of = 0,1 cm for oceanic areas, 
the 10-meter geostrophic wind ratio for neutral conditions 
is, V/Vjo = 0.60* Lettau then allows for stability criteria, 

>vi 

and for the convective case a value of = 0.7^' is arrived 

at. Similarly, =0,5 for the inversion case, 

Advectlon by the frictional wind is derived by first 
decomposing the vector advection as shown by equation (i|.5); 



Then multiplying by O,"^ and rotating the field through the 



W-WTa = + V-Sa 




( 1 * 5 ) 



angle ^ = l5* by a rotational change of coordinates leaves 
equation (46): 
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Assuming the geostrophlc conditions, 66goa = 

T T ® A 

equation (46) may be written in its final form as given by 
equation (47) • 

. 0.9 J. 

•f 

(47) 

^ V / 

Similarly, equation (48) is derived for the advection of 
vapor pressure by t^e 10-meter wind. 

Wlo'WaA = 0.91. Cos lS“X(^o;eA) 

■f 



-0.9lsiV>is°/i«£iSl + 

f I ^ 



( 1 ) 8 ) 



8, Numerical Procedures 

All finite-difference operators used are standard cent- 
ered-diff erences of the t 3 rpe, 



v'^A = laLV^A 






44 ' 






z 



with five-point grids having a mesh distance, d = 38l km 
at 60** latitude, }Iere,l>\ is the map factor for polar stereo- 
graph projections and d is the grid spacing. 



The variables needed for input into the diagnostic omega 
equation were scaled as follows; 
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m 


II 

• 

no 


P 


^ II 

= p. , 2 mb 


y) 


= 9 .2%eo-' 


e 


A q 

= e, .2 mb 


CO 


= Co .2 mb/sec 




A 

= OJJ . 2 cmV sec' 


CJ(.o 


= CcL,,2^ mb/sec 




A ZO 

= 2^.2 cm 


1) 


1\ o'’ 

= Jj .2 cm 




.2 gm-cm*/®6C^ 


T 


= T .2 c 


A 


.4 

= 1 ,J 458 I(. 2 x 10 sin^sec" 



The following physical constants were also needed in the 
computations of Uj , All values of those constants have 
been expressed in cgs units, 

d = 38,1000,000 cm 

^ = 980.0 cra/sec^ 

R = .0287 X 10^ cmVsec^ - “K 

Rv = 0,1461 X 10’ cmVsec^ - *K 

= 1.003 X 10’ cmYsec*- - 

L = 2500 X lo"^ cmVsQC* - ‘'X 

= 0.9771 X lO-*^ “c/cm 

jft> = 0.1600 X 10’*^*C/cm 

^0= 1.213 X 10"^ gm/cm^ 

= 1,055 X 10“^ gm/cm^ 

0.919 X 10'^ gm/cm^ 

The scaled Helmholtz equation, which follows from 
equation (20), expressed in units of, mb - cra*/sec^ , is' 
shown on next page. 
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All height values have been replaced by D-values, the 
deviation of the height field from a standard atmosphere 
height. 

The scaled stability parameter is given by equation 
(50) and is in units of cm*'/sec?‘ 




For computation of the stability parameter D is re- 
scaled to D = D ,2 to prevent overflow when adding the 
standard height, = 290,000 cm. Stability calcu- 
lations are based on 1000-700 mb differences and throughout 
the derivation all gradients and differences at the P| and 
P;^ levels are assumed equal to those at 850 mb and 700 mb 
respective ly. 

To prevent computational instability during the numerical 
Helmholtz solution a minimum value of stability is required. 
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T’lxs valuo copres'aonds to a no.xli’iun lapso of ; 

/ /' A - . V _ 

(To-T, 



V ^ “ 8 






•] 



•2 ' ' 4- 2i ^ooo 
Equation (5l) is the scaled lower boundary condition 
in units of mb/sec. 








< 

Here, the subscript refers to the terrain*height dependent 

-3 

cases, 1, 2, and 3 explained in section (5)» The factor 10 
converts dynes/crn^ to millibars. 



Equations for the diabatic term are computed in terms 
of heating rates, O/Cp , which are assigned to grid points 
according to the existing grid condition: convective, 

inversion, or neutral. 



Equation (52) shows the scaled heating rate for the convective 
case. 






(- 



Pa-?c 



3te>o 




(- 



]AT<oq 









(52) 



'5^300 

The factor, 3600, converts hours to second and heating rate 
is now in units of mb - ‘^C/sec. 



Equation ( 53 ) shows the scaled heating rate for the 
Inversion case (a cooling process) also in mb - **C/sec, 
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The factor X 10^) converts gm-cal/sec to dyne-cm/secj 
K = 24 X 3600 which converts 24-hours to seconds; 0.03 re- 
places 3.0, and permits velocity computation to be made in 
cm/ sec . 



Land area grid points are masked and receive heating 
rate values of zero as do oceanic grid points v;hich satisfy 
the neutral conditions. The remaining oceanic grid points 
receive either heating or cooling values depending on the 
prevailing condition at the point. The diabetic terra is then 
computed and the scaled form in units of mb-cm /sec^ is shown 
by equation (54). 

^ Q - (5h) 

LK Cp 

Pu 

Here, approximates for ease of computation. 

Other equations involved in the computation of the dia- 
betic terra are scaled as follows: 

(1) The air temperature change is units of °C/hr. 

At / V / 

( 55 ) 
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( 56 ) 



(2) The dew-point temperature change in units of ‘*C/hr, 

_ 0.9 ^ ^ CaSIS’T Cio ) €a^ - 

L 

^ il 3 

(3) The wet-bulb temperature change in units of ‘‘c/hr, 

AT^C> _ /NTn I / Z^H<-A/46<3 \/ /KfA _ \ 

At A^ V *■«''’/ / J •? 

A o>^zzU<Sc • 2 

where, )(v.=-0.‘i'>9l RT^'Z V 

(57) 

(4) The condensation level temperature in units of C. 

Tc = Ta - / g 1 V ta -ii) ( 58 ) 

V 0.?(0 ( A / 

( 5 ) The dew-point temperature in units of *C. 

X . t; _ - «a!) (59) 

(6) The condensation level pressure in units of 

• llibars, s^Pci 4- 



mi. 



and 






( 60 ) 



Here T^a, is computed using values of T^ and Tj^ in equation 58. 



Fti = Pa 



ToiZ'’ +z 9 S-<(» 

(ff^*ddk)z<*zoi.l<o 
^ At ' J 
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A ^ A ^ / 

Here is calculated using values of +ATA/At and 

Tj,, + /K Tf>/At in equation,^ ( 58 ) » (59) > (55) » and (56). 



9. The Computer Program 

The Control Data Corporation l60i| digital computer 
t'jras used in this study. It has a core capacity of 32,768 
words of 48 bits each. An operational field of 1,977 
grid points forming a 5l X 47 octagon Inscribed within the 
9*H latitude circle was employed. The grid-mesh is 38 I km 
true at bO**!! latitude. 

Boundary conditions around the octagonal grid were 
determined by the standard subroutines used, Laplacian, 
Jacobian, and all other five-point center-difference 
operations set the grid boundary to zero. The Helmholtz 
relaxation operation set the edge and the next interior 
border point to zero. 



Equation (49l) was solved by a two-dimensional Liebmann 



Here, f\ is the over-relaxation coefficient and the residual 
at any step, Rn , is. 



relaxation technique wherein 
is given 




point 
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The over-relaxation coefficient used was, 

which allowed convergence in approximately 30 scans over 

the 1977 point grid using an initial guess field of zero. The 

Convergence criterion is that the iteration ceases when € 

2 .. 

defined bjr ^ sA -’X fel.TvS bolO'j 3000 r,b-c:'yGec , 

T’'is value corresponds to a vertical velocity of 0,3 X 10 

^ "1. 

rab/sec for a stability of 83 X 10 cm /sec (the standard 
atmosphere stalDility ) . 

Total computation time was approximately 1 minute and 
30 seconds with 23 seconds required for the heating terra 
and 26 seconds for the boundary. condition. 

Due to the abrupt cut-off criteria used for delineating 
tho three cases of diabatic heating (inversion, convective, 
neutral) the final G/c field was smoothed with a five-point 
smoother of the form, 

A = A + kVa 

The smoothing coefficient k is a constant value of (I/ 8 ) 
and the field was smoothed twice. The smoothing operation 
removed small scale irregularities caused by the cut-off 
criteria of the heating term but also reduced peak values, 

A standard PIT;’/P filtering process followed all operations 
involving the Laplacian, .This process removes all small 
scale features with wave numbers greater than 1 $ at latitude 

wh 
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10. Results and Conclusions 

A series of four successive 12-hourly data-sets be- 
ginning with the OOZ maps of 28 April 1966 and ending with 
12Z, 29 April 1966 were used. Since the only diabatic 
heating mechanism introduced into this study was that aris- 
ing by "conduction" from the underlying oceanic surface, 
the results have been depicted only over the North Atlantic 
Ocean, This region has a large density of reporting ships 
and thus the reason for its selection. 

The computational omegas and their associated parameters 
for OOZ 28 April 1966 are contained in figures (6) through 
(11) which are appended. Of these, not only the diabatic 
omegas are shown, but also, the additive effect of the dia- 
batic influence as contrasted with adiabatic computations. 

In figures (12) through (17) only the final-product omegas 
of this study are shown together with the corresponding 
PF/JP analyzed 850 mb D-fields to serve as idontlf iors, 
that is, to indicate qualitative coherency botv;oen the verti- 
cal motions and the associated motion systems. 

In the sequence of figures (6) through (11), it is of 
interest to note the pattern similarity blrtween comparative 
adiabatic vert5^cal motion computations .those, produced here, 
and, for the same times, those; produced!^ PlPiP, which are based 
upon the procedure described by Haltlner et al j^9 The 
magnitude and position of the updraft-centers show a strong 
similarity with those of Plf.TP, however, the model described 
by this paper shows larger downdraft areas extending east- 
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ward from Newfoundland, This apparent discrepancy may be 
attributed to the two differing treatments in 
18 and 19). The^U^at the continent-edge essentially becomes 
the boundary condition for the oceanic computations. 

On the other hand, however, one of tho raajor objectives 
of this study was to test a simplified which uses only 
parameters pertaining to the standard levels of analysis 
(equation 23). In this connection, the greatest time- 
consuming aspect of the PNI-JP co -computation is that of 
deriving and associated parameters at terrain height ro- 

ouiring pressure extrapolation. The general similarity in 
the fields of the two operational procedures which 

have been discussed is evident by an inspection of figures 
(l8) and (19), and this can perhaps be a justification for 
continuing the simpler computations. 

The diabatic effects, which were obtained for the OOZ 
April 28 synoptic time, depict values of h/Cp, so that the 
large positive center southeast of Newfoundland is realistic 
with regard to the northerly flow passing over the Gulf 
Stream, The maximum heating value for the four map periods 
over the Gulf Stream is oyCp = 0,06? mb - *C/secv which 
corresponds to a layer mean heating rate of 0.3 '*C/hr for the 
1000-700 mb layer. The equivalent thickness change for this 
layer is approximately 3 meters/hr. 

For comparison, Petterssen's values [_17^ lor the 
1000-500 mb thickness change during a similar synoptic 
situation in late March show a maximum value of 6 to 8 
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meters/hr in the vicinity of the Gulf Stream. These values 
for sensible and latent heating were computed using fliox 
calculations at the surface based on empiricisms similar 
to those of Laevastu. It should be noted that heating rates 
for this time period in late March are indeed greater than 
late April values. 

Finally, with regard to the diabatic term, note that 
in the area of qualitative verification no significant dia- 
batic cooling values were observed. These areas are limM-' 
ed in extent in winter and springy but could be pf greater 
synoptic consequence in the summer and fall seasons. 

The effect of the diabatic term as it appears on the 
forcing function side of equation (2l|) gives rise to figure 
(11), applicable at approximately 8^0 mb, but actually rep- 
resenting the layer-mean A s could be anticipated from 
the standpoint of heat injection into moving parcels, the 
change (inclusion of the diabatic term) has been such as 
to expand the southern rim of the trough which extends to- 
wards the southwest from Iceland. Similar aspects of coher- 
ency between the resultant co of this study and the trough 
movement and development may be traced out, feature by 
feature as one proceeds through the synoptic sequences. 

For example, the pronounced trough over the Atlantic has, 
by OOZ of 29 April, become oriented North-South just east 
of Greenland. At the same time the updraft cell has taken 
on this same orientation just west of the same trough, so 
that theLP-field gives some confirmation of the dynamic 



processes involved in these map changes. 

The model for vertical motion, as presented by this 
paper, appears to be quite sucessful and future plans Involve 
inclusion of a more realistic frictional term. These more 
realistically computed omegas (extended to 500 mb) are then 
to be used for feedback to yield a prognostic z-field, hour- 
by-hour, The ultimate aim, of course, is to realize a small- 
er R,M,S, error as the motion systems progress over the ocean, 
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Figure 7 
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